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Abstract 

The Yang-Baxter equation for a SU(2) x ?7(l)-symmetric S = \ spin-orbital chain 
was solved using the special computer algorithm developed by the author. The 7 new 
i?-matrices separated on 4 groups are presented. Among the obtained integrable models 
there are special cases related to ID ferromagnet TDAE — C§q, ID superconductors ACgo 
(A=K, Cs, Rb), the quarter filled ladder compound NaA^Os and the model of correlated 
electrons on a chain of Berry phase molecules. 

1 Introduction 

At the beginning of 1970-s Kugel and Khomskii pQ and independently S. Inagaki [2] sug- 
gested two various models to describe magnetic properties of solids with orbital degeneracy in 
electron systems of atoms. Starting from the two-band Hubbard model [3] they have obtained 
low-energy Hamiltonians depending on both spin and pseudospin (orbital) operators. Kugel 
and Khomskii took into account geometry of d-orbitals entailing to asymmetry of hopping 
integrals and obtained a general but realistic Hamiltonian. On the contrary S. Inagaki postu- 
lating the symmetric hopping have distinguished between the Coulomb repulsions on the same 
and different orbitals. It was suggested in pQ- [3] that a non-trivial coupling between spin and 
orbital terms would result to a complex magnetic behavior such as a ferromagnet ism induced 
by orbital ordering. 

While the spin dependence of the Hamiltonians p],[2] has purely 5'L r (2)-invariant Heisen- 
berg form its dependence upon pseudospin is more complicated and in the simplest case possess 
only the U{1) symmetry related to rotations along the z axis in the pseudospin space. 
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After the works [lj,[2] some new applications of 577(2) x [/(l)-invariant spin models were 
suggested in a number of papers pE)-[10]. Ground state and excitations of some SU(2) x U(l)- 
invariant spin chains were studied in [llj,[12j. 

Suggested by L. D. Faddeev and his school Quantum Inverse Scattering Method (Algebraic 
Bethe Ansatze) jl3]-[15j is the most elaborated approach for exact detailed analysis of a one- 
dimensional integrable spin chain. The latter is described by a Hamiltonian, 



N 



H — H n , n+ i, 



n=l 



acting on the finite-dimensional space (<C M )® N (M = 2,3,4,..). The each term i7 nri+1 acts 
nontrivially as a M 2 x M 2 matrix H only on the tensor product of two neighbor spaces 



The keystone of this approaches is the Yang-Baxter equation: 

R 12 (X - n)Rn(\)RiM = i?23(At) J Ri2(A) J R 23 (A - /x), 
with the initial regularity condition: 



(2) 



77(0) = cl. (3) 

Here R(\) is a M 2 x M 2 matrix, 7 is the matrix unity while c is an arbitrary nonzero constant. 
If the Hamiltonian density matrix H relates to R{\) by the following formula: 

dR{\) 



H 



A=0> 



(4) 



then the system (1) is integrable. 

In [16] a new method was suggested for solving the Eqs (2), (3). It is based on the series 
expansion for 7?-matrix: 

R(X) = J2^R {n) X n , (5) 



where, 



7^ 

R( 3 
R^ 
R^ 
R^ 6 



n=0 



77, 
H\ 

H 3 + K, 

77 4 + 2(777<" + KH), 

77 5 + L + 2(KH 2 + H 2 K) + QHKH, 

77 6 + T07 3 + 77 3 7sT + 9(H 2 KH + EKE 2 ) + 

+ 107\" 2 + 3(777 + LE). 



(6) 
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The matrices K and L may be obtained from the following integrability conditions |13j.[16].|17j: 
K 23 — K 12 = [H 12 + H 23 , [if 12, H23]], 

L 23 -L 12 = [Hf 2 + Hl 3 + 3(K 12 + K 23 ),J] + 3{H 12 [J,H 12 ]H 12 + H 23 [J,H 23 ]H 23 ) + 
{H\ 2 H 23 + H 23 Hi 2 )(K 23 — K12) + (K 23 — Ki 2 ){Hi 2 H 23 + H 23 Hi 2 ) — 
2(H 12 (K 23 - K 12 )H 23 + #23(^23 - K 12 )H 12 ), (7) 

where J = [H 12 , H 23 }. 

As the Eq.(2) as the Eqs (6) and (7) are invariant under the following Q e SL(M) action: 

X ^Q- 1 (^Q-'XQ^Q, (8) 

(Q e SL(M), X = H,K, R(X). 

It was shown in [16] a detailed analysis of series expansions applied to quotients of the 
i?-matrix entries gives possibility to guess right the whole i?-matrix. An alternative approach 
for solving the system (2), (3) is given in the recent paper [T9] . 

In the next sections we shall present the obtained by our approach /^-matrices related to 
the general £77(2) x [/(l)-symmetric spin-orbit Hamiltonian, 

H n , n+ i = (s„Sn+i)(ai + a2(rX+i + T-n r n+i) + a3<< +1 + la 6 (< + r^ +1 )) + 

+ + ^n+i) + a 5 r*r* +1 + -a 7 (r* + r* +1 ), (9) 

which may be parameterized by the set of coefficients: S = {ai,a 2 , ...,0,7}. Here in (9) the 
spin and pseudospin operators are expressed from the Pauli matrices: 

s k = ^a k ®I 2 , r k = ^I 2 ®a k , k = x,y,z. (10) 

The following change of coefficients: 

{a 1: a 2 , a 3 , a 4 , as, a^, aj} — > {ai, a 2 , a 3 , a 4 , as, — a^, — a 7 }, (11) 

does not destroy an integrability or change the spectrum of Hamiltonian because it corresponds 
to the transformation (8) with Q = I 2 ® a x . . 

For a chain with even numbers of sites the same is true for the following change of variables: 

{ai, a 2 , a 3 , a 4 , a§, a^, a 7 } — > {a 1; —a 2 , a 3 , — a 4 , as, a^, 07}, (12) 

which corresponds to the graduated version of (8). 

H2n,1n+1 4T2 n i72n,2n+lT2n5 ^2n-l,2n ~^ ^T 2n H 2n -l t2n T 2n . (13) 
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The case, 

S KH = {1 - a, 0,4(1 + a), 0,1 + a, 4,-1} (14) 

corresponds to the Kugel-Khomskii model of ID perovskite pQ. 
The case, 

S 1 = {-(«-/?) + 7 , 2(a + (3), 2(2 7 + (3 - a), ~(3/3 - a), i(3/? + a - 2 7 ), 0, 0}, (15) 

(where unphysical region a ~ ~ 7, but (3 > a and /? > 7) corresponds to the Inagaki's 
model [2] . In [1] it was also applied to organic ID ferromagnet TDAE — Ceo and in [5] to the 
family of ID superconductors ACqo (A = K, Cs, Rb) (with T c > 30K) 
The special case of (15) with, 

a = 1 - 5, (3 = 1 + 5, 7 = 1 - 5 2 , (16) 

(where in physical region < S < 1) was studied in [6]. 
The two cases, 

£NaV 2 o s = { ai8 (3,4a,2f3,4(3-a,0,0}, (17) 
£NaV 2 o s = {_ a? 4 a5 4 tt) 3^3^-4^^ (18) 

correspond to limiting cases of the model describing the quarter-filled ladder compound a' — 
NaV 2 5 [7]. 
The case: 

5" = {1,2, 0,0, 0,0,0}, (19) 

corresponds the effective spin-tube Hamiltonian suggested in [8]. 

Let us also mention the paper [9] where the spin-orbital Hamiltonian was applied to arrays 
of quantum dots. 

In the special S"?7(4)-symmetric point: 

S su & = {1,4,4,1,1,0,0}, (20) 

the model (20) was solved in [18] (-R(A) = r)I + XV). This point corresponds to degenerative 
cases of (15) (a = /3 = 7) (or 5 = 0) Using the transformation (12), (13) we may obtain the 
i?-matrix for the model with 

^ (4) = {1,-4,4,-1,1,0,0}. (21) 
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In the special Sp(4)-symmetric point: 

^ = {1,4,8,2,1,0,0}, 



(22) 



(equivalent to S Sp ^ = {1, -4, 8, -2, 1, 0, 0}) related to 3a = (3 = 7 in (17) the i?-matrix was 
presented in [20] . 

Except (20) and (22) no integrable cases of the Hamiltonian (9) were studied up to now. 
In order to start this process we have solved the system (2), (3) related to (9). The calculations 
were performed by two steps. On the first using the Grobner package of the computer algebra 
system MAPLE 7 we have found 7 new solutions of the system (7). On the second we derived 
the corresponding i?-matrices using the approach suggested in [16], [19J. 

For convenience of representation the obtained i?-matrices are separated on 4 groups. In 
the each one all /^-matrices have similar positions of non-zero entries therefore they may be 
presented in a unique form. 

The author is grateful to P. P. Kulish for statement of the problem and to M. J. Martins 
for helpful discussion. 

Everywhere below e = ±1. 

2 The group 1 



In this group the i?-matrix corresponds to 

={0,4,0,1,0,46,4, 

and has the following form, 



(23) 
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(24) 
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where f± = ±eX + r], g x = A, g 2 = -g 3 = e\. 



3 The group 2 

This group consists of two related subgroups. For the first one the f?-matrix is the following, 



R {2a \X) 



I f+ 
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and there are two solutions. The first one, 

S^ = {l,8e, 4,2^,3,4,-1}, 
corresponds to /+ = /_ = A + 77, g\ — —g 2 = sX. The second, 

S (2a,2) = {i )8 £,4,2£,-l,4,3}, 



(25) 



(26) 



(27) 



corresponds to f± — rj ± A, eg\ — g 2 — X. 

For the second subgroup the corresponding Hamiltonians may be obtained by from (26), (27) 
by the transformation (11) while the R- matrices by transposition with respect to the second 
diagonal. 
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4 The group 3 



In this group .R-matrices have the form, 
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There are two solutions. The first one corresponds to 

£(3,i) = {0,4,0,1,0,0,0}. 



(28) 



(29) 



Here /+ = /_ = sinh(A + rj) for 8 = 2 cosh 77 > 2, /+ = /_ = sin(A + 77) for 8 = 2 cost? < 2 
and f + = /_ = A + rj for 8 = 2. The latter solution is related to the special case of (15) with 
a — f3 and 7 = as well to the special case of (17) with a = 0. 
The second solution corresponds to 



S (3,2) = {o ? 4,0, 1,0,0, 6>}. 



(30) 



Here f± = sinh(?7 ± A) for 8 = 2 cosh rj > 2, /± = sm(r] ± A) for 8 = 2 cos 77 < 2 and f± = rj ± A 
for 9 = 2. In both the cases g = sinhA, h = sinliT? for 8 > 2, g = sinA, h = sin?? for 8 < 2, and 
g = A, h = rj for 8 = 2. 
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5 The group 4 



In this group i?-matrices have the form, 



i? (4) (A) 
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(31) 



where g 2 = fi - / 3 - 

There are two solutions. The first one corresponds to 

S^ = {l,4e, -4, -£,1,0,0}, 



(32) 



where fx = / 2 — sinh(A + rj), eg\ = g 2 = sinhA, / 3 = sinh(A + rj) — sinh A, sinhr] = a/3. For 
e = l this model is the special case of (15) with = 7 = 0. For e = —1 it is SU{2) x SU(2)- 
symmetric (in fact ,Sf/(4)-symmetric [21]). and as it was mentioned in [11] corresponds to the 
four-critical point in the phase diagram of the ferromagnetic SU{2) x [/(l)-symmetric spin- 
orbital model. This model was also suggested in |10j as a model of correlated electrons in a 
lattice of Berry phase molecules. It was also shown in [19] that it is a Temperly-Lieb system 
[22]. 

The second solution corresponds to, 



where fi = f 2 e = 4e 2A - 1, f 3 = 2e A + e 



5(4,2) = {2, 4e, -8, -e,5,0,0}, 

gi = e(l — e~ 2A ), g 2 = 4sinhA. 



(33) 
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